Simulations of interference effects in gated two-dimensional ballistic electron systems 
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We present detailed simulations addressing recent electronic interference experiments, where a 
metallic gate is used to locally modify the Fermi wave-length of the charge carriers. Our numerical 
calculations are based on a solution of the one-particle Schrodinger equation for a realistic model 
of the actual sample geometry, including a Poisson equation based determination of the potential 
due to the gate. The conductance is determined with the multiprobe Landauer-Biittiker formula, 
and in general we find conductance vs. gate voltage characteristics which closely resemble the 
experimental traces. A detailed examination based on quantum mechanical streamlines suggests 
that the simple one-dimensional semiclassical model often used to describe the experiments has only 
a limited range of validity, and that certain 'unexpected' periodicities should not be assigned any 
particular significance, they arise due to the complicated multiple scattering processes occurring in 
certain sample geometries. 
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I. INTRODUCTION 

Recent years have witnessed many experimental and 
theoretical advances addressing the physical properties 
of mesoscopic samples, i.e., structures where the phase 
of the electronic wave function directly affects the mea- 
surable properties Q. A standard way to modify the 
phase of the wave function is to use external magnetic 
fields: electrons traversing the sample along a given path 
will accumulate a phase <f> = (e/h) J A • dl, and thus give 
rise to interference phenomena, such as Aharonov-Bohm 
oscillations. Recently Yacoby et al. g], || demonstrated 
another way of affecting the phase: a biased metallic 
gate, placed above the two-dimensional electron gas, will 
change the electron density (or, equivalently, the local 
Fermi wave-length) underneath it, and thus introduce a 
phase difference between electronic paths that pass un- 
der the gate, and those that do not. In the first ex- 
periment H the amplitude of the interference signal was 
used to extract the energy, or temperature, dependence 
of the dephasing length in a ballistic system. The ex- 
perimental findings allowed a detailed comparison with 
theoretical predictions Q thus underscoring the impor- 
tance of this new technique. The second experiment || 
was the first demonstration of a double-slit interference 
experiment in a solid state system. Both of these ex- 
periments were analyzed with the help of the following 
simple model. Assuming that the charge density is con- 
stant under the gate (but different than elsewhere in the 
sample), Refs. 0, Q find that the phase difference A(f> 
of two representative one-dimensional paths is given by 
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w is the width of the gate, &f 
menta of the unmodified 2DEG and the 2DEG under the 



and k' F are the Fermi mo- 



gate, respectively, and Vdep is the gate voltage required to 
entirely deplete the region under the gate ||.) Indeed, 
the measured conductances have a periodic component 
which essentially scales with the square-root of the gate- 
voltage. In spite of this qualitative agreement, some 
outstanding problems remain. In particular, a Fourier 
analysis of the periodic signal of the double-slit experi- 
ment H contained an unexpected low frequency compo- 
nent, approximately at half of the frequency of the main 
feature. Yacoby et al. S tested a number of plausible 
causes for this behavior (geometric effects, higher elec- 
tronic subbands, and spin-orbit interactions), but con- 
cluded that none of these could satisfactorily explain the 
observations, which thus remained a puzzle. 

The semiclassical picture discussed above is very per- 
suasive, and indeed it has been used in a large number of 
other contexts as well. The geometries studied by Ya- 
coby et al. ||, || are quite complicated, and the possi- 
bility remains that an analysis based on one-dimensional 
straight paths may miss some essential features. To the 
best of our knowledge, these structures have not been 
analyzed in terms of a full solution of the appropriate 
Schrodinger equation, and the purpose of this paper is to 
report such a study, the aim being the resolution of the 
problems encountered in the interpretation of the second 
experiment || . Our work can be seen as a natural exten- 
sion of several recent works reporting detailed solutions 
of the Schrodinger equation for experimentally relevant 
semiconductor nanostructures The resulting wave- 
functions often display a very rich structure and even 
surprising physics, and as an example we mention vor- 
tices around nodal points Q . The new ingredient in our 
work is that, in addition to focusing on two recent exper- 
imental geometries, we include the effect of the phase- 
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modulating gate (PMG) on the potential landscape in 
which the electrons move. Further, we generate stream- 
lines from the probability-current flow; this allows us to 
quantify the role played by the various paths contribut- 
ing to conductance. Our main conclusions are as follows. 
The effects of the gate can be felt in large regions in the 
sample, and, in particular, in the double-slit geometry 
the PMG also affects the slit region. It turns out that 
a description based on a few characteristic paths works 
reasonably well in the geometry of the earlier experiment 
; this is not at all the case for the double slit geome- 
try. In general we find that the resulting conductance 
vs. gate voltage curves are very sensitive to the details 
of the geometry. In particular, the shape of the emitter 
and collector quantum point contacts is found to play an 
important role. For certain parameters the simulated 
conductance curves resemble closely their experimental 
counterparts, yet in other cases, with nominally small 
changes in the parameter values, even the qualitative ap- 
pearance can change drastically. The numerical evidence 
suggests that one should not assign major importance to 
specific features in the Fourier transforms of the periodic 
conductance curves: they may just reflect some details 
of the sample geometry and do not allow a simple semi- 
classical interpretation. 

We should note an important limitation of our numer- 
ical calculations. If we express all energies in terms 
a parameter Eq = h 2 /(2m*d 2 ), where d is the width 
of the injecting electrode, the experiments typically in- 
volve energies of the order of 20000-30000 (we estimate 
d from published electron micrographs). Our numerical 
resources do not allow energies much higher than 5000, 
i.e., one-fourth to one-sixth of the experiments. As a 
consequence, the range of gate- voltages we can study is 
somewhat smaller than what can be achieved experimen- 
tally, but nevertheless we believe that our simulations 
have direct relevance on the reported measurements. 

This paper is organized as follows. In Section II we 
describe the method of calculation, Section III is devoted 
to the analysis of the first experiment, and Section IV 
presents our results for the double-slit geometry. 

II. THE MODEL 

The generic structure considered in this work is shown 
in Fig. [|. It consists of the emitter, collector, phase- 
modulating gate (PMG), and (possibly) the double slit 
(DS). The emitter and collector are modeled as quan- 
tum point contacts. An important role is played by 
the four base contacts: electrons scattering off from the 
DS or PMG, and not making it to the collector leave 
the device via these contacts and do not contribute to 
further interference patterns. In the simulations these 
ideal base contacts are represented by open (Neumann) 
boundary conditions. In effect, then, we are consider- 
ing a six-terminal geometry (emitter, collector, and four 
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FIG. 1. Schematic representation of the six-terminal struc- 
ture. The boundaries of the structure and the double slit 
are shown as shaded areas; they are modeled as hard walls. 
The lithographic placement for the phase-modulating metallic 
gate (PMG) which is deposited on the surface of the sample 
a height h above the 2DEG (typically h ~ lOOA) is indicated 
as a white rectangular area. The corresponding potential is 
shown in Fig. 2. 

base contacts). The boundaries defining the structure 
are taken as hard walls, and the potential describing the 
PMG is described below. Our numerical method for the 
solution of the Schrodinger equation is quite standard, 
and here we only give those special features that are nec- 
essary for understanding the computed results, given in 
several figures below. The incident wave function in the 
n-th transverse channel is given by 

ipin,n (x, y) = sin (imx) e~ tk " v 

+ r nm sin (irmx) e lkmV (1) 

m 

where the sine- functions are the transverse eigenfunctions 
of the injecting electrode of unit width {d — 1), r nm are 
the corresponding reflection coefficients, and the wave- 
vectors k n are defined via 

e = kl + 7T V (2) 

where the energy e is given in units of Eq defined above, 
and k n is in units of 1/d. Analogously, the collector 
wave-function is expressed as 

Vw,n {x, y) = ^ tnm sin (irmx) e ~ lkmV , (3) 
m 

where t nm is the transmission coefficient from mode n to 
mode to, and the collector is assumed to have the same 
width as the emitter. Finally, the boundary conditions 
at the s-th base contact (of width L s ) are specified by 

V> s ,„ (x,y) =^t 5 , nm -^sin L m V -^l\ e <*..«(*-*.). 
m V s \ s / 

(4) 
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FIG. 2. The characteristic shape of the electrostatic poten- 
tial due to the phase-modulating gate. Also are shown curved 
emitter and collector contacts as used in the experiment of 
Ref.[2]. 
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FIG. 3. Conductance vs. gate voltage for the device of 
Fig. 2 (shown as inset) for initial energies e = 2500 (top) 
and e = 5000 (bottom). The predictions of the semiclassical 
formula for conductance maxima are shown as arrows. 



Here x s , y s are the coordinates of the walls defining the 
base contact s, and kg m — e — ir 2 m 2 /L 2 . 

The Hamiltonian in a tight-binding representation for 
the two-dimensional electron system is 



H 



-*(4,j c i.i+i + c l,i c i+i,J + h - c -) + V id4,j l 



(5) 



where Vy is the electrostatic potential due to the PMG. 
We use the transfer matrix method as formulated in Ref. 
H to compute the various transmission and reflection 
coefficients (and hence the conductances) ||. 

A compact theory for the electrostatic potential caused 
by metallic gates on the surface of the heterostructure 
has been developed by Davies et al. Jl(| for a number 
of different physical assumptions (pinned surface, frozen 
surface, linear screening etc.). Following their analysis 
we choose |llj (here h is the distance between the gate 
and the 2DEG): 



V, 



V(x,y,z = h) = [f (x, y - w/2, h) 
-f(x,y + w/2,h)], 



where 



f (x,y,h) — arctan 



R — x — y 



(6) 



(7) 



R = \/x 2 +y 2 



h 2 . 



A typical potential profile is shown in Fig. 



III. THE EFFECT OF PHASE-MODULATING 
GATE 

We first consider the experiment reported in The 
aim there was to study dephasing due to electron-electron 
collisions in a ballistic sample, and the phase modulating 



gate was introduced to generate an interference signal: 
the amplitude of the oscillatory component of the con- 
ductance is a direct measure of the phase-coherent part 
of the electrons, and thus allows one to extract the phase- 
breaking rate as a function of injection energy or temper- 
ature, and compare it to theoretical predictions The 
actual sample had two 'semi-infinite' gates, however only 
one of them was activated and we therefore model this 
system with the geometry shown in Fig. |^. Let us first 
construct an analytic estimate for the expected behavior. 
Fixing the coordinates so that the PMG runs parallel to 
the x-axis, and that the electron moves along the y-axis, 
the semiclassical formula for the phase accumulated un- 
der the PMG is 



9 = i J dy^2m* (E -V (x,y,h)). 



(8) 



Neglecting the paths that pass near the edge of the PMG 
we can take the x — > oo limit of Eq.(0), and find 



V(y,h) = ^ 

7T 



arctan 



y + w/2 



arctan 



w/2 - y 



which, to a good approximation, provided that h <^ w, 
can be approximated by a rectangular barrier of width w 
and height V g . The phase difference between electrons 
passing under the PMG and those which do not is then 
readily found to be 



A6 = w(^- V /T 



(9) 



(Again, all energies are in units of Eq.) This derivation 
provides a justification for the phenomenological expres- 
sion for phase A0 used in the experimental articles [0, 
||, see also Jl2]. Maximal constructive interference oc- 
curs when A9 = 2nn, which leads to expected maxima 
for PMG biases at 



(n) = 4ir\/en/w — (2irn/w) z 



(10) 
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FIG. 4. The modulus of the wave-function for two gate 
voltages V g = 0, 1250 for incident electron in the first subband 
of the emitter with energy e = 2500. 

Fig. [}] shows the numerically computed conductance 
vs. PMG bias for incident electron energies e = 2500 
(corresponding to 16 propagating modes), and e = 5000 
(24 modes). The conductance was evaluated with the 
Landaucr-Buttikcr multichannel formula lOl, |14|: 



2e 2e \ - . |2 k m 

nm 



(ii) 



The voltages at the base contacts are set to zero. The ar- 
rows in Fig. H indicate the predictions of the semiclassical 
formula, and one observes a good overall agreement. As 
could be expected, the conductance does not show any in- 
terference structure when the gate voltage becomes very 
large, because the effective potential due to the gate is 
then nontransmitting. 

In order to gain a deeper understanding of why the 
one-dimensional model works so well in this particular 
case it is useful to study the quantum mechanical stream- 
lines (for applications to several other physical systems 
one may consult, e.g., @, @, @, |J], @). We use 
the following construction. Writing the wave-function in 
terms of an amplitude and a phase, 



t/j — -y/pexp (iS/h) 



(12) 



the real and imaginary parts of the time-independent 
Schrodinger equation, 



yield 



2m 
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FIG. 5. Stream-lines for the gate voltages of Fig. 4. The 
stream-lines passing under the gate and ending in the col- 
lector are essentially straight lines and can be described by 
one-dimensional quantum mechanics. 
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(16) 



where 



According to Bohm one interprets the electrons as 
"real" particles in the classical sense, following a contin- 
uous and causally defined trajectory with a well-defined 
position x with the momentum given by m x= V5. The 
force acting on the particle is not derivable from the 
classical potential V alone, but acquires a quantum me- 
chanical contribution from Vqm, Eq.(16). The current 
stream-lines can then computed as in classical mechanics 
but including the quantum force. The stream-lines can 
be viewed as an alternative graphical presentation of the 
quantum mechanical probability current density, see, e.g. 
Fig. 1 of Ref. @. 

Figure |] shows the computed wave-function for two 
different values of the gate potential V g , the right panel 
corresponding to the first maximum in Fig. ||, top panel. 
We direct attention to the following features, (i) The 
wave-functions display a rather regular pattern even at 
a finite gate voltage (which breaks the mirror symmetry 
of the problem), (ii) The curvilinear injector leads to 
a clear focusing effect. (iii) The effective wave-length 
is clearly longer under the gate than elsewhere in the 
2DEG, in accordance with the expectations. Fig. || 
shows the computed stream-lines. We note that most of 
the stream-lines ending in the collector are, even in the 
case of a finite gate potential, essentially straight lines. 
Thus they can be described by one-dimensional quantum 
mechanics, and consequently the semiclassical estimate 
for occurrence of conductance maxima, shown as arrows 
in Fig. ^, works reasonably accurately. We next turn to 
the double-slit geometry, where matters turn out to be 



4 




50D 1000 1500 2000 2500 3000 3500 4000 4500 5000 



FIG. 6. Conductance vs. gate voltage for double-slit geom- 
etry, with curved emitter and collector. The top view of the 
device is shown as an inset. 
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FIG. 7. Conductance vs. gate voltage for double-slit ge- 
ometry, with rectangular emitter and collector. A schematic 
view of the device is shown as an inset. 

quite different. 

IV. DOUBLE SLIT GEOMETRY 

We now introduce the double slit. In order to compare 
most directly with the results obtained in the previous 
section, we first consider same emitters and collectors as 
before, even though the experiment was done with a dif- 
ferent design (this will be discussed below). Fig. || shows 
the computed conductance in the presence of the double 
slit, as indicated by the inset in the top panel. Again, 
we see oscillations in the conductance, however the val- 
ues of the gate voltage at which the conductance is at 
maximum do not correspond to the values predicted by 
the simple estimates, such as Eq.(H). We next consider 
the experimental geometry of Ref . B , where the emitter 
and collector are rectilinear. Figure [?| shows the com- 
puted conductance curves and Fig. ^| displays examples 
of computed wave-functions. Conductance oscillations 
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FIG. 8. The wave-function for double-slit geometry at two 
gate-voltages, V g = and V g = 2250, the latter correspond- 
ing to the last significant maximum in the conductance vs. 
gate-voltage characteristic of top panel of Fig. 7. 

are quite evident in Fig. [?], but it is much harder to find 
any regular periodicity, in contrast to the curves for the 
device without the double slit of Fig. [| Since our simu- 
lated conductance curves only have few maxima (because 
of the computational restrictions to relatively low ener- 
gies) we did not find a Fourier-analysis helpful (as was 
the case for experiments which had a larger available gate 
potential range) : the resulting spectrum is dominated by 
spurious edge effects. Finally, Figures ^| and [II] show 
the computed stream-lines for the curved and rectan- 
gular emitters/collectors, respectively. We draw atten- 
tion to the qualitatively different picture as compared 
to the device without the double-slit: the approximately 
straight-line form of the stream-lines is almost entirely 
lost. Most importantly, the stream-lines passing under 
the gate show a rather irregular structure with a very 
wide range of effective path-lengths under the gate. It is 
instructive to consider the pair of paths denoted by A and 
B (which have symmetric initial conditions at the emit- 
ter) : the combined effect of the double-slit and the gate 
is to distort B significantly, and it is not surprising that 
an one-dimensional model fails to describe it properly. 

The device with rectangular emitter/collector, as in 
S, has another intriguing property: there is a clear ten- 
dency to form resonances between the emitter and the 
double-slit (and, to a lesser extent, between the double- 
slit and the collector). The effect of these resonances 
can be understood in terms of a beating phenomenon: 
their frequency mixes with that due to the PMG, and in 
general one can expect much more irregular conductance 
vs. gate voltage curves, as is the case with curved emit- 
ters/collectors. It is quite conceivable that this mixing 
can contribute to the half-frequency oscillation observed 
in pi. Another indication of these 'size-resonances' is that 
the conductance is not always maximum at zero gate volt- 
age (see, e.g., the lower panel of Fig. 0): this is because 
the emitter quantum point contact is not always matched 
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to the resonator modes of the cavity formed by the gates 
defining the emitter and the double-slit, and a finite gate- 
voltage can move the resonator modes so as to achieve 
more efficient injection from the emitter. In view of our 
simulations it would appear to be interesting to repeat 
the double-slit experiment with curved emitter/collector: 
the experimental trace is expected to be easier to inter- 
pret because one achieves a better focused injection and 
diminishes complications due to the resonator modes. 

V. CONCLUSION 

We have presented simulations of phase-coherent 
charge transport in gated mesoscopic structures. The 
simulations can describe the experiments at least quali- 
tatively, and under certain circumstances quantitatively. 
Using quantum stream-lines as an interpretative tool 
we are able to offer an explanation of why certain ex- 
periments can be interpreted with the help of one- 
dimensional models, while others cannot. We find in gen- 
eral great sensitivity to geometric effects, however these 
can be controlled at least to some extent by careful de- 
vice design aided with simulations of the kind presented 
here, in particular when extended to higher energies. 
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FIG. 10. Stream-lines for double-slit geometry with rect- 
angular emitter and collector. 
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